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ABSTRACT 

The hard state of X-ray binaries (XRBs) is characterized by a power law spectrum in the X-ray band, and a 
fiat/inverted radio/IR spectrum associated with occasionally imaged compact jets. It has generally been thought 
that the hard X-rays result from Compton upscattering of thermal accretion disk photons by a hot, coronal 
plasma whose properties are inferred via spectral fitting. Interestingly, these properties — especially those from 
certain magnetized corona models — are very similar to the derived plasma conditions at the jet footpoints. 
Here we explore the question of whether the 'corona' and 'jet base' are in fact related, starting by testing the 
strongest premise that they are synonymous. In such models, the radio through the soft X-rays are dominated 
by synchrotron emission, while the hard X-rays are dominated by inverse Compton at the jet base - with both 
disk and synchrotron photons acting as seed photons. The conditions at the jet base fix the conditions along the 
rest of the jet, thus creating a direct link between the X-ray and radio emission. We also add to this model a 
simple iron line and convolve the spectrum with neutral reflection. After forward-folding the predicted spectra 
through the detector response functions, we compare the results to simultaneous radio/X-ray data obtained 
from the hard states of the Galactic XRBs GX 339-4 and Cygnus X-l. Results from simple Compton corona 
model fits are also presented for comparison. We demonstrate that the jet model fits are statistically as good as 
the single-component corona model X-ray fits, yet are also able to address the simultaneous radio data. 
Subject headings: X-rays: binaries — black hole physics — radiation mechanisms: non-thermal — accretion, ac- 
cretion disks — X-rays: general 
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1. INTRODUCTION 

Bipolar plasma outflows, loosely termed jets, are a feature 
common to a variety of different astrophysical objects. Jets 
are observed in accreting compact objects of all scales, as 
well as at stellar birth and death. They seem equally able 
to form out of massive accretion disks as well as from quasi- 
spherically collapsing plasma, which suggests that their cre- 
ation is a basic byproduct of some routinely occurring ingredi- 
ents. These necessary inputs seem to be rotation, collapsing or 
infalling/accreting plasma, and magnetic fields. Jets likely are 
formed at least in part from a nec essity to shed the system o f 
excess angular momentum (e.g., Meier & Nakamura 2004). 
The exact details of their formation, as well as the nature of 
the coupling between the infalling/collapsing plasma and the 
outflow is still an active area of controversy and research on 
many fronts. However, the fact that so many types of systems 
share this mechanism extends hope that clues gained from one 
class of objects can be applied to the others in the search for a 
unified picture. 

The advantage of looking for answers in X-ray binary 
(XRB) jets is twofold. First, their formation seems to be a 
recurrent phenomenon on time scales we can observe repeat- 
edly over the course of our own human lifetimes. For in- 
stance in the Galactic XRB GX 339-4 and Cyg X-l, accre- 
tion without observable jets proceeds into accr etion with dis- 
cernible jets and back a gain every few year s ( Fender ll2003t 
iPottsch midt et al] l2003t iGleifiner et aTJ l2004at iHoman et alJ 
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2005). Thus the same object can provide tests of theories de- 
veloped based on its earlier activity. Secondly, the existence 
of quality simultaneous broadband data for many sources 
linking the "traditional" jet outflow bands (radio and likely 
also infrared, or IR) and inflow bands (X-ray and optical) 
gives us a way of studying the relationship and interplay be- 
tween the two sides of accretion. We use the word "tradi- 
tional" to refer to the picture which existed up until recently 
that the jets and accretion flow radiate in different energy 
bands and can generally be studied as distinct phenomena. 
This picture has rapidly evolved in the last few years, begin- 
ning with the discovery of an unexpectedly intimate connec- 
tion between the low and high frequen cy wavebands in the 
hard state (HS) of XRBs (see, e.g. jMcClintock & Remillard! 
2003, for a detailed definition). 

Using radio observations in conjunction with 1.3-12.2keV 
data from the All Sky Monitor ( ASM) on-board th e 
Rossi X-ray Timing Exp lorer (RXTE; Levine et al. 1996), 
Hannikain en et all {1998) first discovered a near-linear cor- 
relation between the radio and X-ray emission of GX 339-4. 
These studies were followed by joint rad io observations and 
pointed X-ray obse rvations with RXTE ( Wilms *et al J 1 1 9991 
Nowak et al. 2002) . This intensive simultaneous monitor- 
ing showed tha t, in fact, the corre lation is non-linear, with 
L R oc L^ 7 ( Cor bel et afll2000l 120031) . Additionally, this non- 
linear correlation holds over orders of magnitude changes in 
the source luminosity with time. In fact, even after fading into 
quiescence, GX 339-4 seemin gly returns to the s ame correla- 
tion in subsequent outbursts (Corb elet alJ. 12003). although it 
is not clear whether the correlat ion always ma intains the same 
normalization in each instance ( Nowa k et alJ2005l) . 
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Despite some possible variation in its normalization, the 
correlation between radio and X-ray fluxes appears to be 
universal to all H S XRBs with comparable broadband data 
JGallo etalJ 120031) . The 0.7 slope of the log(L R )-log(L x ) 
correlation follows directly fro m analytic predictions of 
scaling synchrotron jet models ( Falcke & Biermann 119951 
iMarkoff etai]l2003l) . and can be generalized for other X-ray 
emission p rocesses in terms of the ir dependence on the accre- 
tion rate (Heinz & Sunvaev 2003). The correlation's normal- 
ization in the (Lx,Lr) plane also depends on the central engine 
mass. When this scaling is accounted for, unbeamed, extra- 
galactic, supermassive black holes sources agree remarkably 
well with the same radio/ X-ray correlation found i n Galactic, 
stellar-mass sources (Merloni et al. 20031 iFalcke et al.12 004). 
The underlying physics governing the connection between the 
radio/IR and X-ray bands must therefore be fundamental to 
accr eting sources regardles s of mass. 

In Markoff et al. (2001a) we showed that jet synchrotron 
emission could account for the broad continuum features 
of the simultaneous radio through X-ray hard state obser- 
vation of the Galactic XRB, XTE Jl 118+480. In a later 
work, we showed that this same model could also explain 
the broad spectral features of thirteen simultaneous or quasi- 
simultaneous radio/X-ray observations of GX 339-4, and that 
the radio/X-ray correlation was natural ly produced by only 
changing the power input into the jet (Markoff et al. 2003). 
These models, however, did not attempt to address the fine 
features present in the X-ray spectrum which are the hall- 
marks of an optically thick accretion disk: in particular fluo- 
rescent iron lines and a characteristic flattening above lOkeV 
attributed to Compto n reflection Cl.ightman & White! 119881 
iRoss & FabianlfT993l) . The reflection is generally assumed 
to result from a hard X-ray continuum originating above the 
cooler accretion disk. 

We began to explore the interacti on of the jet em ission with 
the putative accretion disk in Markoff & Nowak (2004). In 
the hard state (HS), the spectrum is hard enough and the 
observed reflection signature weak enough to be problem- 
atic for coronal models where there is significa nt coverage 
of the disk by hard X-ray-emitting material (e.g., iDove et all 
and references therein). Various mechanisms have been 
proposed to decrease both the cooling of the hard X-ray- 
emitting material by soft disk photons and the fraction of 
refle cted X-rays. T hese mechanisms include patch y coro- 
nae (ISter n et all! 995l). high disk ionizati on ( IRoss et al J 19991 
Navakshin 2000; Ballant vne et al.ll2001l) and beaming of the 
coronae away from the disk wi th mildly relativistic velocities 
(Belo borodovll999llMalzac et alJ20 01). This latter approach 
seems extremely close in principle to the characteristics of the 
base of a jet. 

Compact accretion disk coronae are theoretical concepts 
based on observations that suggest the presence of hot elec- 
trons radiating near the inner parts of a thin accretion disk. 
The existence of hot electrons in the same region is also em- 
pirically required by the jets we image in the radio wavebands. 
The high brightness te mperatures an d occasionally measured 
linear polarization (e.g., Fender 2001 ) argue strongly that non- 
thermal synchrotron is the dominant radio-emitting process. 
Because the base of the jet is significantly more compact than 
the outer region which dominates the cm-bands, conservation 
arguments imply an even hotter, denser medium near the ac- 
cretion disk. We therefore wish to determine whether the base 
of the jets can "subsume" t he role of the corona . This idea was 
proposed already back in iFender et al.l (119991) . and by sub- 



sume we mean provide the spectral characteristics (both direct 
and reflected) which are traditionally attributed to a compact 
Comptonizing corona. Thus we will compare coronal fits us- 
ing a standard Comptonizing model (see below) to fits made 
by our jet model. Regardless of the ultimate relationship de- 
termined to exist between the base of the jets and the corona, 
this is a critical step in the road to understanding. The results 
of this study will provide valuable clues as to the nature of 
the corona, both in terms of geometrical as well as internal 
characteristics which differ from the base of the jets. 

As an additional note, the concept of cyclo-synchrotron 
photons feeding inverse Co mpton processes in a corona has 
existed for over 20 years ( Fabian ^t alJ Il982h . inspired by 
multiwavelengt h observation s showing similar optical/X-ray 
behavior ( Motch et al] ll982|). It has recently been reconsid- 
ered for several sources l|di Matteo et alJl!997llMerloni et al.l 
2000; Wardzihski & Zdziarski 2000) in the context of mag- 
netic flares in the corona. While a valuable step towards un- 
derstanding the role of magnetic fields, these models do not 
address the radio emission from the jets with the same model. 

The remainder of this paper is structured as follows. In 
Section [2] we very briefly summarize the model, with a full 
description included in Appendix lAl In Section[3]we present 
our data analysis and fitting techniques, the results of which 
will be presented in Section^ We discuss our conclusions in 
Section|5] 

2. MODEL DETAILS 

2.1. Jet Model Direct Emission 

Our model for the jets in HS XRBs is based upon five main 
assumptions: 1) the total power in the jet scales with the to- 
tal accretion power at the inner edge of the disk, Mc 2 , 2) the 
jet is expanding freely and is only weakly accelerated via the 
resulting pressure gradient, 3) the jet contains cold protons 
which carry most of the kinetic energy while leptons do most 
of the radiating, 4) particles are eventually accelerated into a 
power-law distribution, and 5) this power-law is maintained 
along the length of the jet thereafter. These assumptions are 
motivated by observations which are discussed further in Ap- 
pendix |A| where we have also included a detailed descrip- 
tion of the model and its development. In this section we 
only summarize the key parameters and assumptions of the 
jet model. 

Beyond a small nozzle region of constant radius, which es- 
tablishes the base of the jet, the jet expands sideways with 
the initial proper sound speed for a relativistic electron/proton 
plasma, 7 s /3 s c ~ 0.4c. The plasma is weakly accelerated along 
the resulting longitudinal pressure gradient, which allows an 
exact solu tion for the v elocity profile via the Euler equation 
(see, e.g., [Falcke| 119961) . This results in a roughly logarith- 
mic dependence of velocity upon distance z. After a period of 
more rapid acceleration immediately beyond the nozzle, the 
velocity gradient lessens, and the velocity saturates at Lorentz 
factors of Tj >2-3. The size of the base of the jet, ro, is a free 
parameter and once fixed determines the radius as a function 
of distance along the jet, r(z). 

The the model is most sensitive to the fitted parameter Nj, 
which acts as a normalization. It dictates the power initially 
given to the particles and magnetic field at the base of the jet, 
and is expressed in terms of a fraction of the Eddington lumi- 
nosity LecM = 1.25 x 10 38 Mbh.o erg s" 1 . The total power input 
at the base of the jets is in fact approximately an order of mag- 
nitude larger than A^Ledd, due to the requirement of distributed 
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acceleration along the jets (see Appendix). Once Nj is spec- 
ified and conservation is assumed, the macroscopic physical 
parameters along the jet are determined. We assume that the 
jet power is evenly shared between the internal and external 
pressures, and that the radiating particles are close to equipar- 
tition with the magnetic field. The radiating particles enter 
the base of the jet where the bulk velocities are lowest, with 
a quasi-thermal distribution. Around 10 - 100 r g , a significant 
fraction of the particles are accelerated into a power-law tail 
whose maximum energy is determined self-consistently with 
the local cooling rates. 

The particles in the jet radiatively cool via adiabatic expan- 
sion, the synchrotron process, and inverse Compton upscat- 
tering; however, adiabatic expansion is assumed to dominate 
the observed effects of cooling because of distributed accel- 
eration. While thermal photons from the accretion disk are 
included as seed photons in our Compton calculations, the 
beaming reduces their energy density compared to the rest 
frame synchrotron photons (synchrotron self-Compton; SSC), 
except at the very base of the jet where they can be of the 
same order. Reprocessed disk radiation will contribute even 
less and thus its feedback on the X-ray spectrum will be neg- 
ligible. We therefore do not include this latter component in 
our calculations. In the case of near-equipartition of energy 
between the magnetic field and particles, inverse Compton 
processes only will dominate the direct synchrotron emission 
close to the compact object where relativistic beaming is at its 
minimum. 

2.2. Accretion Disk Feedback Signatures 

As the jet plasma is only weakly beamed, some fraction of 
the radiation will impact upon the cooler material in the accre- 
tion disk an d lead to fluorescen t line emiss ion and a reflection 
hump (e.g.. iLi ghtman & Whitdfl988h . In Marko ff & Nowakl 
( 2004) we calculated the resulting reflection spectrum for typ- 
ical jet models, and found that when synchrotron emission 
from the acceleration region dominates the hard X-ray spec- 
trum, the fraction of reflected emission is very small (~few 
percent). In contrast, when Compton processes (predomi- 
nantly SSC) originating near the base dominate the hard X- 
rays, 'reflection fr actions' of < 20% a r e poss ible. 

As discussed by Markoff & Nowak (2004), calculating the 
reflection spectrum for our jet models is complex, as one 
needs to calculate a reflected spectrum associated with each 
emission region along the axis of the jet. Furthermore, for 
our detailed jet models the disk sees a qualitatively differ- 
ent jet spectrum than the direct spectrum viewed by a distant 
observer. This i s in contrast to the simple 'moving corona' 
models of, e.g., iBeloborodovl (l999), where both the direct 
spectrum and the spectrum impinging upon the disk have the 
same simple power law shape, and the total reflected spectrum 
can be calculated analytically. For the initial spectral fit stud- 
ies presented here, a detailed reflection calculation for our jet 
model was deemed too computationally prohibitive. 

For the purposes of the study presented here, we adopt a 
simplified approach. In our fits, we add a single Gaussian 
line to our continuum model, and we allow the line energy to 
vary between 6-7 keV and the line width, a, to vary between 
0-1.5keV. We then convolve the entire directly viewed spec- 
trum with a non-rela tivistic reflection model derived from the 
Green's functions of Magdziarz & Zdziarski ( 1995). The am- 
plitude of the reflected spectrum is left as a fit parameter, ex- 
pressed in the usual manner as a fractional solid angle, £1 /2ir, 
subtended by the reflector. (Q/2tt= 1 indicates that the ra- 



diation impinging upon the reflector is equal to that directly 
viewed by the distant observer.) This general approach of phe- 
nomenologically fitting a reflection spectrum is similar to that 
employed by many p ure Comptonization models (e.g. lCoppH 
1 1 999t iPoutanenl 1 9981) . albeit here it is less likely that the di- 
rectly viewed spectrum is a completely adequate proxy for the 
spectrum viewed by the cold disk. We will consider more so- 
phisticated approaches to calculating the reflection features in 
future works. 

2.3. Fitting Method 

In our previous papers (e.g.. iMarkoff et aLll2001al 12003). 
we only compared the jet model to spectra that had been 'un- 
folded' with XSPEC using simple exponentially cutoff, bro- 
ken power-law models. This is in fact the usual practice 
among many researchers when comparing complicated theo- 
retical models to mul tiwavelength data (e.g., Esin et alll997t 
Mark offet al.ll200Tah . Unfortunately, such a procedure does 
not allow the determination of a statistical goodness-of-fit, 
and more importantly it does not allow one to compare the 
model to the fine features of the X-ray spectral density distri- 
bution (SED). Here we address this issue by forward-folding 
the jet model through the detector response matrices of the 
X-ray instrument, and then comparing to the data in 'detector 
space' . 

Specifically, we have imported our jet + multi-color black- 
body thermal disk (sim ilar to diskbb; Mitsuda et al.lll984c 
Makishima et all 19 86) model as a subroutine for use in stan- 
dard X-ray data analysis packa ges. The c urrent routine works 
in bot h XSPEC JArnaudll9 96) and I S I S llHouck & Denicolal 
2000). In this work, we perform our fits using I S I S v. 1 .2.6 as 
this latter package can read in and fit lower frequency simul- 
taneous data sets (e.g., radio through optical data) as ASCII 
files, without the need for creating dummy response matri- 
ces. Furthermore, the ISIS 'unfolded spectra' (shown in the 
figures throughout this work) are independent of the assumed 
spectral model (i.e., the unfolding is done solely with the re- 
sponse matrix and effective area files; see iNowak et alf2 005 
for details). All shown residuals, however, are for compar- 
isons between the pro perly forward- folded model and the data 
in detector space. See Nowak et al. (2005) for a discussion of 
further differences between ISIS and XSPEC. 

We are making detailed comparisons between our jet model 
and detector space data for the first time. Indeed, to our 
knowledge, this is the first time that such a complex, multi- 
wavelength model has been compared properly to the X-ray 
data in detector space. For this paper, we therefore have cho- 
sen to initially leave most jet model parameters free so that we 
can fully explore their effect on the fits. We can then deter- 
mine which parameters are most fundamental for describing 
the data and which can effectively be frozen in future appli- 
cations. This is essentially in the same spirit as prior studies 
with Compton ization models, where models such as eqpair 
( Coppil IT999I) have significantly more parameters than are 
used in a typical fit. For such Comptonization models, even 
fairly basic parameters, such as thermal vs. nonthermal com- 
pactness (i.e., the ratio of coronal energy to coronal radius), 
can be degenerate with one another in real data fits, and one 
often chooses to freeze the nonthermal parameters to negligi- 
ble values. For a detailed disc ussion of these points, see the 
review article by Coppi ( 2004|). 

Currently, our multi-component model runs significantly 
slower than single-component models which calculate Comp- 
tonization in the corona. Determining meaningful error bars 
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Table 1. Log of Cyg X-l and GX 339-4 RXTE Observations 



Source 


ObsID 


Date 


Exposure 


PCUoff 








(sec.) 




Cyg X-l 


40099-01-19 


1999 September 25 


8080 


13 


Cyg X-l 


60090-01-26 


2003 February 23 


9696 


1/4 


Cyg X-l 


60090-01-41 


2003 September 22 


11600 


13/14 


GX 339-4 


20181-01-02 


1997 February 10 


10528 




GX 339-4 


40108-02-01 


1999 April 2 


9152 




GX 339-4 


40108-02-03 


1999 May 14 


9776 





NOTE. — The table lists those proportional counter units (PCUs) which 
were off during (part of) the observation. A slash denotes a logical "or". 



for all parameters was not always possible because for some 
fit parameters the resulting spectrum can vary noncontinu- 
ously. Given that this is the first time that broadband data are 
being properly compared to our forward-folded model, such 
initial difficulties are not surprising. Nevertheless, the model 
fits presented here represent an exploration of a significant 
amount of the parameter space, especially given the fitting 
time scales involved. 

3. DATA 

We model the data from Rossi X-ray Timing Explorer 
(RXTE) observations for the hard state Galactic black hole 
candidates (BHCs) GX 339-4 and Cygnus X-l. We use 
data from both the RXTE Proportional Counter Array (PCA; 
Jahod a et alJU996|) an d from the High Energ y X-ray Timing 
Experiment ("HEXTE: iRothschild et aT1ll998l) . The data have 
been extracted with the recent HE AS OFT 5.3.1 software 
release. Compared to earlier releases of this software, the rel- 
ative calibration of the PCA and the HEX TE is now in ex - 
cellent agreement (see analysis in, e.g.. IWilms et al.l 2005). 
Power-law fits to spectra of the Crab pulsar and nebula show 
the PCA data to have a systematic uncertainty of 0.5%, which 
we added in quadrature to the data. We used the energy bands 
3-22 ke V for the PCA (and grouped the data to a minimum 
of 30 counts/bin) and 1 8-200 keV for the HEXTE (where we 
added data from the two clusters and grouped the counts to 
a minimum signal-to-noise, after background subtraction, of 
10 in each bin). We use the top layer of the PCA only. Dur- 
ing part of the observations, some proportional counter units 
of the PCA were switched off. In these cases we added the 
spectra extracted for the different PCU combinations, appro- 
priately weighting the response matrices. 

The Cygnus X-l data consist of three of the many observa- 
tions taken during our several year s long RXTE monitoring 
campaign (Pottsch midt et alJ 120031 iGleifiner etaD I2004albl 
and references therein) and were chosen to be representative 
spectra bracketing the typical hard state spectral variations of 
Cyg X- 1 . The d ata from GX 339-4 also have been discussed 
by us elsewhere (IWilms et alll999tlNowak et al 120021120051) . 
and bracket the typical variations in the overall hardness and 
luminosity of the GX 339-4 hard state. (They do not, how- 
ever, encompass the most extreme bright and relatively soft 
HS data discussed bvlHoman et alJl2005lNowak et aljl2005l 
and B elloni et alJl2005l ) Table Ulcontains a log of the obser- 
vations. 



In order to compare our jet models to the more standard 
Comptonization models, we also fit our observations with the 



hybrid thermal/non-thermal Comptonization model eqpair 
(Coppil 119991) . We slightly changed the publicly available 
version of this code to make the model run under XSPEC 
v . 1 1 . 3 . 1 , and hence also under ISIS v. 1.2. 6. To 
allow for a partial covering of the disk by the Compton 
corona, we add emission fr om a disk using the diskpn 
model JGierliriski et aljfl9* 99). The peak temperature of this 
disk model is also set equal to the peak temperature of the 
seed photons input to the Compton corona. These seed pho- 
tons are similarly assumed to have a spectral energy distribu- 
tion that follows that of the diskpn model. The corona is 
presumed to be comprised solely of thermal electrons, and its 
properties are described by two parameters: a seed electron 
scattering optical depth, t p (pair production can increase the 
total optical depth from this value, although pair production is 
negligible for the fits described in this paper), and the coronal 
compactness (i.e., coronal energy content divided by coronal 
radius), expressed as a ratio relative to the soft, seed photon 
compactness, £f,/£ s - 

We also allow for an additive Gaussian Fe Ka emission 
line whose centroid energy is constrained to be between 6 
and 7 keV and whose width is constrained to be a < 1.5 ke V. 
In our use of the eqpair model, the reflection compo- 
nent is relativistically smeared using the velocity profile of 
the diskpn disk model, which in turn follows the rota- 
tional velocity of a standard thin disk circularly rotating un- 
der the influence of a Paczyhski-W ita pseudo-Newtonian po- 
tential (Pa czvhskv & Wiitafl9~8 0). The reflecting medium is 
assumed to have solar abundances; however, the reflecting 
medium can be ionized, as pa rameterized by the ionization 
paramet er, | (Don e et alll992l) . 

Table lA4l gives the best fit parameters of these eqpair 
fits, and the fits are shown in the third column of Figs.|3]&El 
The reduced x 2 values are close to unity, indicating very good 
agreement between the model and the data. The parameters 
found for the Comptonization model are in general agreement 
with those reported for earlier observations. The fits yield 
moderate optical electron depths of r es ~ 1.5-3 with coronal 
compactnesses of Ih/ts ~ 3.5-8, which correspond to coronal 
temperatures of ~ 25-100 keV. (Note that the r es ^6.7forObs 
ID 40108-02-03 has rather large error bars, and is partly being 
driven by the poor statistics in the HEXTE band.) 

All of the coronal model fits yield reflection covering frac- 
tions < 20%. None of the fits strongly require ionized re- 
flection; however, the reflection model used by the eqpair 
code is very sim plified (being based upon the pexriv code; 
iDone et alJU992l) . The fitted equivalent widths of the Fe Ka 
lines are on the order of 150 eV, and consistent with the highly 
smeared reflection models employed, the lines are generally 
found to be rather broad. As will be discussed elsewhere in 
more detail (Wilms et al. 2005), we believe that part (although 
not all) of the fitted line broadness is caused by the Comp- 
tonization model employed here having too simplified a de- 
scription of the transition between the accretion disk emission 
and the Comptonization continuum. Finally, the changes in 
accretion disk temperature and flux, coupled with the changes 
in coronal compactness and optical depth, are in response to 
the fact that we have deliberately chosen observations that 
span a wide range of flux and spectral hardnesses. 

Detailed discussions of coronal model fits to Cyg X-l data 
(e.g., comparison to phenomenological broke n power-law 
fits, co mparison to timing data) are discussed in Wil ms etlri] 
(2005). Their are three main points that we wish to em- 
phasize with the fits presented here. First, the eqpair fits 
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describe the data well - as well as any coronal models that 
we have explored (although not better than simp le exponen- 
tially cutoff, brok en power-law models; iNowak et alJl2005l 
Wilm s~e7al] 120051) . Second, the fits do not require compar- 
atively as large reflection fractions 2 . Third, and related to 
the previous point, the coronal models themselves describe 
a continuum spectrum more complicated than a simple, ex- 
ponentially cutoff power-law. Some fraction of the spectral 
hardening above lOkeV, normally attributed to the 'reflection 
hump', is in fact described by curvature of the Comptoniza- 
tion continuum. In Compton coronal models, this curvature is 
in part due to transiting from a spectrum strongly affected by 
the soft seed photons to one strongly dominated by photons 
that have been multiply Compton scattered. As we discuss 
below, this continuum curvature an additional interpretation 
in jet models. 

4. JET MODEL RESULTS 

As mentioned above, our strategy in testing our jet models 
against the data was to allow the maximum number of model 
parameters to vary freely during the fit. This allowed us to 
explore which parameters had the greatest influence on the 
fits. Several parameters settled on fairly similar values for 
both sources and all observations, which suggests that those 
parameters could be fixed in future applications of this model. 
We will discuss this possibility explicitly below. 

Fits with a complicated model can easily lead to false, local 
minima; therefore, we began the fitting procedure outside of 
ISIS, using unfolded data sets in order to find a set of starting 
parameters that would yield reduced x 2 ^ 5. For this paper we 
have chosen to explore only models with rough equipartition 
between the magnetic and radiating particle pressures (k « 1). 
This assumption is obviously not applicable for Poynting flux 
dominated jets; however, we consider the very weak accelera- 
tion required by observations to be an indication that magnetic 
domination (which would generally imply stronger accelera- 
tion mechanisms) is not likely to be very extreme . 

We use the results of Markoff & Nowak (2004) as a rough 
guide for consistency, and only consider models with reflec- 
tion fractions < 20% for jet nozzle radii ro < 20r g to be 'suc- 
cessful' fits. For smaller values of the radii, slightly lower 
values for the reflection fraction would be expected. In all 
cases the jet models seemed to naturally prefer reflection frac- 
tions n/2ir < 20%, so in practice this restriction on reflection 
fraction was not an issue. 

Fig. [2 shows a representative model for GX 339-4 in or- 
der to illustrate the contribution from the various components 
which go into the later figures showing actual fits. The ra- 
dio through IR originates exclusively from self-absorbed syn- 
chrotron radiation beginning at z acc and continuing outwards 
along the jet. Optically thin synchrotron emission from the 
accelerated power-law tail of leptons also contributes to the 
soft X-ray band. The base of the jet radiates direct syn- 
chrotron photons, giving a slight hump in the optical/UV, 
which are then included as seed photons (along with disk 
photons) for upscattering by the emitting electrons into an 
SSC/EC "hump" in the hard X-rays. The shape of this hump 

2 The pexr iv model, incorporated into the eqpai r model, tends to have 
an unnaturally sharp ionized Fe edge structure, esp e cially when compared 
to the more sophisticated models of Ross & Fabian 12005). The reflection 
fraction can be artificially depressed when employing the pexriv model, 
so as to minimize residuals from its sharp ionization edge. However, as we 
employ a large degree of relativistic smearing, and hence smooth out the 
edge, we do not expect to be subject to this systematic effect. 



is rounded, due to the quasi-thermal particle distribution as- 
sumed in the base, and thus reduces the need for a large 
fraction of disk-reflected photons to contribute to the spectral 
break/hardening at ~ 1 ke V. The possibility that there are two 
correlated continuum components in the X-ray band is sup- 
ported by our ensemble of simultaneous radio/X-ray observa- 
tions of Cyg X-l, whose spectra are ver y well-described by 
exponentially cu t-off, broken power-laws jNowak et al.l2005l 
Wilms et al1 l2005l) . In these models, the soft X-ray spectral 
slope is very well-correlated with the break of the hard X-ray 
spectral slope, and the hard X-rays are very well-corr elated 
with the radio (IdleiBner et alJ2004atlNowak et alJ2005l) . The 
multicolor blackbody included in the jet model is also shown 
in the soft X-ray band. The photons from this component are 
included in the inverse Compton upscattering within the jet. 

The statistical fits are shown in Figs.|2]&[3] for GX 339-4 
and Cyg X-l, respectively. Each figure encompasses nine 
panels, with each row representing a single simultaneous ra- 
dio and X-ray observation of the source. The first column 
shows the entire radio through X-ray jet model + soft disk + 
reflection fit, and the second column focuses just on the X-ray 
band. The last column shows the thermal Compton corona 
model from eqpai r for comparison. Values for all fitted pa- 
rameters are given in Tables 2-4. 

The results presented in Figs. |2] & E and Tables IA2I IA3I 
& IA4l immediatelv demonstrate three important results of this 
work. 

• Jet models describe the data equally well as pure 
Compton corona models, even when employing 
the broad-band and high statistics of the RXTE 
data. There has bee n question in the literature (e.g., 
Zdzi arski et al J 120041) as to whether or not jet models 
can adequately describe the spectral cutoff of the hard 
tail. While pure synchrotron jet models may not be 
able to describe the steep cutoff present in some (but 
not all) observations, jets can also successfully account 
for this hard tail cutoff via Comptonization. The pri- 
mary differences between 'traditional' coronal models 
and jet coronal models are that the coronal tempera- 
tures are higher in the latter (~ 10 I0 K vs. ~ 10 9 K), 
and the disk seed photons in jet models are signifi- 
cantly augmented — if not dominated — by synchrotron 
seed photons. 

• The jet models, similar to the Compton corona mod- 
els, describe a complex continuum with curvature 
- specifically, a hardening above 10 keV. Again, the 
hardening often predominantly ascribed to a 'reflec- 
tion hump' is being somewhat subsumed by contin- 
uum emission. In the jet model, the steeper slope of 
the soft X-rays is primarily due to synchrotron emis- 
sion from the jet, rather than being influenced by 
the disk photons, while the hardening is due to the 
SSC/EC component. Reflection is present, but rep- 
resents a smaller fraction of the total hard flux com- 
pared to the corona models. Note that the fits dis- 
cussed here are in contrast to earlier incarnations of 
the jet model which attempted to describe observations 
solely via the synchrotro n component (e.g., th e stud- 
ies of XTE Jl 118+480; Mai-kof f^tail 12001 ah . For 
XTE Jl 118+480, the lack of any discernible break or 
hardening near lOkeV in the continuum spectrum was 
used to argue for the lack of a reflection component 
(iMilleretaU 12002). which would be consistent with 
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FIG. 1. — Representative model for GX 339—4 obs.id. 40108-02-01, showing jet and multi-color blackbody models only. This is the model which is then 
convolved with the disk reflection and line emission elements for the final statistical fit. The components are labeled. 



a jet spectrum dominated b y the synchrotron compo- 
nent jMarkoff et alJl2001at Markoff & Nowakl [2004). 
In contrast, the fact that all of the XRB spectra shown 
here harden above lOkeV argue for the importance of 
both the synchrotron and SSC/EC components of the 
jet spectrum for these fits, as well as the presence of 
reflection. 

• Although the jet and Compton coronal models de- 
scribe the X-ray data equally well, only the jet model 
naturally describes both the amplitude and slope of 
the radio data without the need for additional free 
parameters. It is important to note that the radio spec- 
trum, as well as the bulk of the 3-200 keV X-ray con- 
tinuum (i.e., the overall amplitude and continuum cur- 
vature), are almost solely driven by the parameters of 
the jet: energy input to the jet (iVjLEddX radius of the jet 
base (ro), electron temperature at the jet base (T e ), slope 
of the power-law component of the electron energy dis- 
tribution (p), the equipartition parameter (k), and the lo- 
cation of the particle acceleration zone in the jet (z ace ). 
If one begins instead with a corona model and attempts 
to fit the amplitude and slope of the radio data by as- 
suming a relationship to an extant jet model, this will in 
essence require the addition of many more parameters 
to describe just two physical quantities. We argue that 
nothing can be gained by this approach, which is why 
we have chosen to study the test case of full coupling. 

From our set of jet model fits that successfully describe the 



GX 339-4 and Cyg X-l data, we can now explore the impli- 
cations of the fitted parameters. Cyg X-l and GX 339-4 are 
somewhat different sources while in their hard states. The lu- 
minosity of GX 339-4 varies much more than that of Cyg X- 
1 . This is apparent in the fitted values for the power normal- 
ization N\ and output power Lj shown in the first two columns 
of Table IA2I Both sources show best fit Nj in the range of 
10" 4 - 10" 3 L Edd , but GX 339-4 shows significantly more vari- 
ation in this parameter. Differences can also be seen in the 
other main parameters, for instance on average GX 339-4 fa- 
vors slightly larger values of the nozzle radius, ro. This re- 
flects the slightly higher ratio of X-ray to radio flux observed 
in Cyg X-l. A smaller scale jet base, for a given jet power, 
increases the jet compactness which in turn gives a slightly 
higher flux (with the X-ray flux being more sensitive to this 
effect) and pushes all jet emission to slightly higher frequen- 
cies. 

Other fit parameters also indicate differences between 
these two sources; however, the parameters surprisingly fall 
roughly in the same range given the a priori possibility for 
much greater parameter variations. Note that the lowest lumi- 
nosity GX 339-4 observation, Obs ID 40108-02-03, has poor 
enough statistics that the following general statements about 
the GX 339-4 fits will be based mostly on the brighter two 
observations. For example, both the electron temperature, T e , 
and the location of the start of the acceleration zone, z acc , tend 
to be larger in GX 339-4 than in Cyg X- 1 . (Although we were 
unable to find adequate error bars for the latter parameter, the 
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trend of finding larger values in GX 339-4 was persistent.) 
Additionally, we note that the electron power law index, p, 
shows more variation in GX 339-4 than in Cyg X-l. 

Although the range in electron temperatures at the base of 
the jets falls within a range of a factor of < 2 for both sources, 
one of the criticisms of this class of models has been that ~ 3- 
5 x 10 10 K is not as "natural" a value as the ~ 100 keV typi- 
cally used in thermal Comptonization models. The jet model 
electron temperature, however, is comparable to or slightly 
greater than the value typ ically derived for ra diatively ineffi- 
cient accretion flows (e.g., Narava n et al.ll998l) . If some frac- 
tion of the accreting plasma is fed directly into the jets, and 
also perhaps heated slightly in the process, we would expect 
such a temperature. 

On the other hand, several jet model free parameters 
seemed to settle quite quickly into similar values for both Cyg 
X-l and GX 339-4: the equipartition factor k, the fraction of 
accelerated particles plf, and the ratio of nozzle length to ra- 
dius ho. Physically, this suggests that these jets either are close 
to equipartition, or we have found a local minima of solutions 
for k « 1-2, and we may find other reasonable minima for 
k ;§> 1, e.g., magnetically-dominated jets. We plan to explore 
this particular question in a future work. It seems reason- 
able, however, to freeze the fraction of accelerated particles at 
~ 75% and make the statement that reasonably efficient accel- 
eration is expected to occur in hard state BHC jets. Similarly, 
a compact base/corona region with scale height similar to the 
radius seems a reasonable assumption. It is likely that this 
parameter can also be frozen in future applications of the jet 
model. 

We listed the two "acceleration parameters", the shock 
speed relative to the bulk plasma flow, u. dcc /c, and the ratio of 
the scattering mean free path to the gyroradius, f sc , as separate 
fit parameters. However, as discussed in the Appendix, these 
parameters are perhaps not physically meaningful as currently 
defined, since we are no longer as convinced that the accelera- 
tion process is diffusive Fermi acceleration. These parameters 
enter into the acceleration rate as the factor / = (« a cc/c) 2 / f sc , 
which is compared to the sum of the cooling rates from syn- 
chrotron and inverse Compton radiation, and adiabatic expan- 
sion in order to calculate the local maximum accelerated lep- 
ton energy. Therefore, this parameter can be loosely inter- 
preted as a factor related to the efficiency of the acceleration 
process and merged for future fits into a single parameter. For 
the current fits, the individual components were left free to 
vary, with mixed results. For Cyg X-l,/w2- 3x 10~ 4 , sug- 
gesting a meaningful range of fits. However, for GX 339-4, / 
ranges over two orders of magnitude, and we suspect that we 
did not fully explore parameter space meaningfully for this 
source. In future applications of the model, it is likely that 
we will combine these two parameters into a single parameter 
that will absorb our uncertainty about the acceleration pro- 
cess. 

While we list best fit values for the multicolor blackbody 
disk model parameters Ldisk and as well as the derived 
inner radius r m , it is important to keep in mind that because the 
data only extend down to a few keV, we cannot meaningfully 
constrain this component from spectral fitting. As described 
in the Appendix, the disk photons are not as important to the 
overall photon field as the locally produced synchrotron pho- 
tons and thus can mainly be constrained by their direct spec- 
tral contribution. This contribution, although weak, is in fact 
require d for a good fit, but it is not unique. Ibragimov et al. 
(2005) have also fit some of the X-ray data presented here, and 



for similar reasons fixed the thermal disk parameters. On the 
other hand, for the disk emission plus Comptonization models 
prese nted here (as well as for those presented in Wil ms etafl 
2005), the temperature of the disk was tied to the temperature 
of the seed photons — i.e., the thermal photons were the only 
source of seed photons for Comptonization. Thus by virtue of 
these imposed restrictions, much more stringent formal limits 
for the disk components were attainable. 

The main result to take away about the accretion flow mod- 
eling is that our fits are generally consistent with a sub- 
Eddington accretion disk with temperatures somewhat less 
than 1 keV. Similarly, the total power entering into the jets 
is roughly consistent with being of the same order as the ob- 
served luminosity required in the disk to be consistent with the 
data. The presence of a weak disk component is necessary for 
a good fit, and thus disk photons will contribute to the inverse 
Compton component from the jet. We do not, however, feel 
we can confidently make any statements about the disk geom- 
etry and thus the accretion rate at the inner radius assumedly 
feeding the jets. 

For both GX 339-4 and Cyg X-l, the amount of reflection 
required by the jet models is roughly comparable to, although 
in general slightly lower than, that required by Comptoniza- 
tion models. This trend can be understood by the fact that 
the jet base SSC/EC component dominates above 10 keV, and 
to some extent subsumes the role played by the "Compton 
hump" normally attributed to disk reflection. This effect of 
course brings up an interesting point: there is a clear degen- 
eracy in how the spectral hardening above 10 keV can be un- 
derstood in terms of continuum models. The fact that the jet 
SSC/EC component has a similar appearance to the Comp- 
ton reflection hump does not preclude the presence of both. 
Clearly, the fluorescent Fe line implies that there must be a 
degree of reflection. What these results (as well as the Comp- 
tonization fits; Wilms et al. 2005) do suggest is that one can- 
not uniquely determine a reflection fraction independently of 
the presumed continuum model. 

5. DISCUSSION AND CONCLUSIONS 

The key result of this work is that we have clearly demon- 
strated that even the high statistics of RXTE X-ray data cannot 
distinguish between the jet model and thermal Comptoniza- 
tion. (Indeed, the RXTE data also cannot distinguish between 
thermal and nonthermal Comptonization models, as has been 
previously discussed by, e.g., Coppi 2004). In this sense, the 
base of the jets can be said to effectively "subsume", at least 
spectrally, the role created for the corona. The main difference 
between these two pictures then comes down to geometry, in- 
flow versus outflow, and the relationship of the corona to the 
lower-frequency-radiating regions of the outer jets. The jet 
model presented here is consistent with the phys ical picture 
seen in relativistic MHD simu lations ( Miller & Stone] l2000t 
lHawlev & KrolikfeOOlllMeier & Nakamural2004l) . where the 
corona is not static but instead is a windy hot material blow- 
ing away from the inner regions of the accretion flow. This 
wind/corona is likely what enters into, or is collimated and 
becomes, the jets and therefore it is intimately related to the 
larger scale outflows, similar to what we propose here. In ad- 
dition, the jet base provides a natural illuminator for the disk. 

Our ability to constrain the role of the accretion disk is sig- 
nificantly limited by the lack of low-energy data. Ideally we 
would like to be able to determine the inner temperature, and 
perhaps also assess the need for a recessed inner edge, explic- 
itly from the observations. Having this information would al- 
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low us to better constrain the ratio of external to synchrotron- 
produced photons in the jet base. Additionally, we expect 
differences between jet models and Compton corona models 
to become more pronounced as one considers data above the 
nominal Comptonization cutoff at 1 00-200 keV. From both 
of these perspectives, the observational situation will likely 
soon improve with the launch of ASTRO-E2, which will have 
the ability to observe from 0.5-600 keV and thus more strin- 
gently constrain both jet and Compton corona models. Ad- 
ditionally, although the broad-band effects of reflection, e.g., 
the reflection hump, cannot be fit independently from the as- 
sumed continuum, there is some hope that finer spectral reso- 
lution (perhaps coupled to correlated studies of the continuum 
variability) will be able to more uniquely to determine the re- 
flection fraction based upon detailed modeling of the Fe line 
JYoung & Revnoldsl2000tlRos s & Fabi anl2005l) . 

In the scenario considered here, the nozzle/base of the jet 
is the dominant region for the creation of the hard X-ray 
emission via Comptonization processes. Constructing a com- 
pletely realistic description of this region, however, is diffi- 
cult given the that the physics of jet formation, acceleration, 
and collimation are not yet understood. The nozzle essen- 
tially represents the initial conditions which fix parameters 
along the rest of the jet. We have shown that the form we 
have developed for this region over the last several years can 
successfully mimic the corona in purpose, but it is not as well- 
determined feature as we would like. An alternative, but also 
currently rather ad-hoc, possibility would be to invoke spe- 
cial conditions linking an outflowing corona with more typ- 
ical thermal plasma conditions to an outer jet. This would 
require some kind of acceleration and/or heating of the parti- 
cles at the interface (and thus more free parameters), but may 
satisfy many of the spectral requirements. 

One possible way to discern between these two situations, 
or at the very least limit the possible flux from a nozzle, could 
be to calculate the amount of local ionization expected from 
nozzle UV photons and compare that to what is observed (J. 
Krolik, priv. comm.). If it can be shown for a given source that 
all observed line emission can be accounted for by the known 
companion, this limit may be rather stringent. We note, how- 
ever, that in High Mass X-ray Binaries such as Cyg X-l, 
the companion luminosity, especially in the UV, is extremely 
high. It is therefore more likely that meaningful limits can be 
obtained for Low Mass X-ray Binaries, such as GX 339-4. 

While we have focused on modeling simultaneous radio/X- 
ray data sets in this paper, for GX 339-4 and several other 
sources observations do exist which in addition have simul- 



taneous IR/optical data. Near-IR data in particular can pro- 
vide a valuable constraint on the jet synchrotron component, 
especially when the transition turnover fro m optically thick 
to thin is resolved (e.g. lHoman et alJl2005l) . Using the work 
presented here as groundwork, we will explore elsewhere 
whether more parameters can be constrained by such data 
sets. Similarly, with HESS already online and GLAST im- 
minent, predictions of 7-ray fluxes will be another important 
way of testing not only the current model, but models which 
incorporate hadronic processes. 

However, the next major frontier for application of our jet 
model clearly is to consider questions of time-dependence. 
By studying the lag between the radio and X-ray correlated 
variability, we can constrain the plasma speed and/or dis- 
tance between the two regions. The existing radio/X-ray flux 
correlations constrain the long time scale, while simultane- 
ous radio/RXTE poi nted observations c onstrain thousand of 
second time scales llGleiBner et alJl2004al) . We are currently 
using long duration (several hundred ksec), quasi-continuous 
radio/X-ray observations of Cyg X-l to explore intermediate 
time scales (Markoff et al., in prep.). The models considered 
here, however, are steady-state and focus on the hard accretion 
state only. Some of the most interesting and revealing behav- 
ior is seen during transitions between these states. For in- 
stance, during transition into the hard state we observe a hard 
X-ray tail before the jet forms as a det ectable radio structure 
JNowak et aUl2002l [Miller et aT1 l2004l) . The presence of the 
tail suggests that the jet base/corona may form as a viable re- 
gion before the outer stable radio-emitting structures are built 
up. There are also significant clues comi ng from noise and 
faster scaled time variability such as QPOs Homa n~& Bellonil 
(2005) which need to be integrated into the picture. In general 
the state transitions reveal many interesting properties which 
are clearly not due to a steady state model, but which do share 
several similar properties to the more persistent hard states. 
Studying the time-dependent behavior of these systems will 
likely be an effective way to further constrain the conditions 
at the jet base/corona. 
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APPENDIX 

DETAILED DESCRIPTION OF THE JET MODEL 

Background 

Like all models for optically thick outflows, our jet model builds on the initial work of lBlandford & KoniglHl979l) . Their paper 
demonstrated how a superposition of idealized, self-absorbed, conical jet components results in a flat synchrotron spectrum in 
the radio bands. In reality, compact jets usually deviate from this idealized case, and show a slight spectral inversion in the radio 
wavebands, with spectral index a ~ 0.0-0.2 (F v oc v a ). A more realist ic treatment of the internal physics can ac count for this 
inversion. Using the framework for hydrodynamical jets developed in Falck e & Biermannl d 19951) . iFalckd ( 1199 6) developed a 
"free jet" model which is more self-consistent. For instance, the assumption that the jet is roughly conical via lateral expansion 
with a constant sound speed requires weak longitudinal acceleration. The Euler equation can then be solved for the longitudinal 
velocity gradient, assuming only adiabatic losses, and is analogous to a pressure driven wind solution with the proper jet speed 
replacing the wind speed. This weak acceleration combined with relativistic beaming effects will result in the spectral inversion 
typical of compact jet sources. 

These earlier models were very useful for studying the general physical scalings predicted by compact jet models. However, 
without full particle distributions and radiative transfer, they cannot be tested against observed spectral data. Our models have 
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FlG. 2. — Jet and corona model fits to radio and X-ray data of GX339— 4, with individual spectral components for the jet model shown in grey. In each panel, 
synchrotron from the jet pre-shock region is shown as a dot-dash line, synchrotron from the post-shock region is shown as a dash-triple dot line, synchrotron 
self-Compton plus external Compton is shown as a solid line, and the disk radiation is shown as a dotted line. The panels are laid out as follows. From bottom to 
top (left to right on the figure), the fits are with the mf f jet model (radio/X-ray and X-ray data) and then the eqpair model. From left to right (top to bottom 
on the figure), the Obslds are 20181-01-02, 40108-02-01, and 40108-02-03. The spectra are shown as model-independent unfolded spectra; however, the fits and 
residuals are from the properly forward-folded models. 



been developed on the basis of these earlier papers, but specifically focus on the problem of spectral predictions. These models 
are only appropriate for systems where the jets are not expected to be highly collimated nor highly accelerated, e.g., the jets in 
hard state XRBs and LLAGN. 

Jet Physical Parameters 

We start by assuming that the base of the jets are at the speed of sound for a relativistic gas with adiabatic index 4/3, giving a 
proper sound speed /3 S 7 S ~ 0.4. Furthermore, we assume the jets expands laterally with the constant sound speed, which is appro- 
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FIG. 3.— Same as Fig.0 except for Cyg X-l data. From left to right, (top to bottom on the figure), the Obslds are 4099-01-19, 60090-01-26, and 60090-01-41. 



priate under the assumption of maximally efficient jets (see Falcke &Biermann 1995). The resulting longitudinal acceleration 
along the axes is largest near the bases, and asymptotically approaches a maximum of Lorenz factor 7j ~ 2 - 3 in the outer radio 
emitting re gions, and y, < 1-5 for the synchr otron emission. Such weak beaming factors are suggested by observations of hard 
state jets (Maccarone 20031 iGallo et al]l2003l) . in contrast to the near light-speed velocities implied by superluminal motion in 
higher luminosity radio-emitting states. 

The main free parameters of this model are the power normalization Nj , the radius of the jet nozzle ro, the equipartition of energy 
between the radiating leptons and the magnetic field k = Ub/ U e , and the parameters determining the initial lepton distribution 
as described below. We also left the length of the nozzle region (with fixed radius), ho free, but it turns out the fits are not very 
sensitive to its value. Other parameters which come in, such as the central mass my,, the inclination of the jets to the line of 
sight 9\ and the distance to the source d^, are determined from observations and remain fixed for any given source. Finally 
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there are two parameters related to the presence of weak disk emission, the total luminosity radiated by the disk and its inner 
radius temperature, Ld and 7i n respectively. Thes e cannot be well con strained by our data because although they are spectrally fit 
assuming a multi-color disk blackbody (Mit suda et alJll984llMakishima et alll986l) . the low-energy X-ray spectrum is not well 
covered in our observations. Similarly, although they do contribute to the photon field upscattered by energetic jet leptons, they 
are not the dominate component. 

The most important free parameter is iVj, parameterized in terms of the Eddington luminosity, which determines the power 
initially input into the particles and magnetic field at the base of the jets. The total power input at the base of the jets is in 
fact approximately an order of magnitude larger than AfjLEdd, and for our model can be estimated only a posteriori from the 
resulting emission in the emitting rest frame. The difference between N^L^ and the total output radiation Lj expresses the 
lack of knowledge of what initially collimates the n ozzle, as well as w hat causes the distributed reenergization — and consequent 
energy requirement — to maintain the observed (e.g. lJester et alJl2001l) particle distributions along the jet (see below). Nj plays 
a similar role to the compactness parameter in thermal Comptonization models, which parameterizes the total energetics in the 
absence of an understanding of the mechanism for energizing the corona (iDove et alll997tlCoppill999l) . 

As a consequence of assuming maximal jets, there should be approximate equipartition between the kinetic and internal energy. 
Once iVj is fixed, therefore, the energetics along the jets are fixed as well. For a given ro and k at the base, conservation of energy 
and particle flux together with the proper velocity profile along the jet, 7j(z)/3j(z), determine the radius, density, and magnetic 
field profiles along the entire jets. What is left to be determined is how the energy in the radiating particles is distributed. 

Radiating Particle Distributions 

As mentioned earlier, we assume that the jets contain both leptons and hadrons, but that the hadrons remain non-relativistic 
and serve only to car ry the bulk kinetic energy. The leptons, which are assumed to be mildly relativistic and quasi-thermally 
distributed (see, e.g., Ouataert & Gruzinov 1999) are thus more likely to be accelerated. In any case, we assume that leptons 
are predominantly responsible for the observed radiation. This is begging the question of hadronic acceleration, because the end 
effect of significant hadronic acceleration and subsequent inelastic collisions will invariably be relativistic leptons. The expected 
energy distributions, whether from direct acceleration or creation via hadronic collisions, are not trivial to discern from X-ray 
emission alone. Therefore we do not specify whether the radiating particles are electrons or positrons, since this will not effect 
the observable outcome. Eventuall y the consistency with these results can be c hecked against hadronic models, but this is a very 
difficult problem (see, e.g. lMannheimll993tlMarkoff et alJll999tlBottcher & Reimerl200 4) and is not considered here. The one 
hint for what will be required beyond neutral matter is that we do not require n e /n p = 1, and once a fit is made, we can see if 
additional pair creation is necessary. 

The current model was developed in simplified form by Falcke & Markoff (2000), where we first incorporated simple particle 
distributions for the radiating leptons and calculated the radiative transfer along the jets. Our goal was to determine if the same 
kind of model which could explain the inverted radio core of the Galactic Center supermassive black hole, Sgr A*, also predicted 
significant X -ray emission. This idea was motivated by the first identification of Sgr A* in the X-ray band by the Chandra 
Observatory (Baganoff et al. 2003). We explored two "canonic al" types of particle distributions, power-laws and Maxwellians, 
for Sgr A* in its quiescent state (see, e.g., Melia & Falcke 2001 ). These two ostensibly different distributions result in similar fits 
as long as the characteristic par ticle energ y, (7 e , m i n for the power-law or 7e, pe ak for the Maxwellian) is similar. Unlike "typical" 
low-luminosity AGN (LLAGN; Hol l 19991) . the quiescent Sgr A* does not show any indication of a power-law of optically thin 
synchrotron emission after the break from its flat/inverted radio spectrum. Therefore if the radiating particles have a power-law 
distribution, it must be so steep as to be indistinguishable from a Maxwellian in the optically thin regime. Either scenario strongly 
suggests, for example, that acceleration in the jets of Sgr A* is absent or very inefficient at L ~ 10~ 9 LEdd- 

Based on Sgr A*, we therefore assume that the baseline particle distribution, below a certain accretion rate which must be 
very low, is quasi-thermal. This scenario is consistent with particles being advected into the jets directly from the accretion 
flow, or created and roughly thermalized near the base. Sources with jets at higher luminosity, such as AGN and XRBs, do, 
however, show the standard indicat or of particle acceleration: optica lly thin synchrotron emission from a power-law distribution 
of particles (e.g., Marscher & Ge ari l 1 9 85t [Fen der & Kuulkers 2001, and refs. th erein). I nterestin gly, Sgr A* shows d aily flaring 
in which the X-ray spectrum rapidly increases in flux and hardens dramatically (Baganoff et al. 2001; Baganoff 2003). The flare 
emission stems from energization of the particles via acceleration, heating or both at once (Markoff et al. 2001b), resulting in 
enhanced synchrotron and/or SSC emission i n the X-rays. B ecause the increased X-ray flux during the flares brings Sgr A* more 
in line with typical levels of LLAGN activity (Markoff 2005), one interpretation is that acceleration of particles is the main factor 
contributing to the comparative weakness of Sgr A*. We also therefore assume that above some critical accretion rate above that 
of Sgr A*, particle acceleratio n will be commonplace in the jets. This is also necessary to explain the high synchrotron efficiency 
in XRBs (see. e.g.. lMarkoff et al.l2003l) . 

In the only XRB wher e the op tically thick-to-thin turnover from compact jet synchrotron has been directly observed (GX 339-4 
; Corbel & Fen derl2002HHoman et al 120051) . it occurs at a significantly lower frequency than would be expected if the accelera- 
tion began in the jet base. We therefore assume that when particle acceleration occurs, it begins at some location z acc along the jet, 
which is a free parameter but which appears to fall in the range of a few 10— 10 2 r g for most sources we have so far considered. 
The fraction of particles accelerated into the tail is a free parameter in our model, but is generally driven to fairly high values so 
we will likely fix it at some value above 80% in future applications. The power-law index of the tail p {N{E) oc E~ p ) is also a free 
parameter, and can be constrained by the data. 

Particle acceleration must compete with radiative and adiabatic cooling to energize the quasi-thermal particles into a power- 
law tail. Where the sum of the cooling rates equals the acceleration rate defines the maximum achievable particle energy, 7 e ,max- 
Unfortunately since th e ac celeration process is sti ll an open question, the acceleration rate is open to significant interpretation. In 
iMarkoff et alJ (1200 1 bh and Marko ff" st alJ (1200 lal) . we used acceleration rates appropriate for diffusive shock acceleration in the 
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most conservative case of the magnetic field direction parallel to the shock normal (see Joki piill987l) . The rate of diffusive shock 
acceleration is always proportional to the magnetic field. If synchrotron cooling dominates the cooling term, 7 e , m ax will thus 
be independent of the magne t ic field and can give information about the plasma parameters if the location of the cutoff can be 
determined. In Markoff et al. (2001a) we showed that this results in a synchrotron emission cutoff around 100 keV for (3^ > 0.5 
and if the ratio between the particle's mean free path for diffusive scattering to the gyroradius, £ = A/Vb « 100. In general £ 
is thought to be limited by the ratio of the particle to shock velocities, and so likely no higher than 10 2 - 10 3 , but this value is 
certainly not fixed. 

At the time of writing, other processes such as stochastic resonant acceleration (e.g., Miller 1998) are coming into favor. Given 
this uncertainty in the mechanism, we are loosening the constraints on the acceleration rate we used previously and are letting £ 
and f3 s h vary. If the acceleration rate we are using is not really physical, then these parameters have lost their meaning and become 
essentially fudge factors absorbing the free parameters controlling the acceleration rate for other mechanisms, to be determined 
later. If shock acceleration does still hold, then these tell us about the plasma conditions in the jets. 

Regardless of the process, the inferred cooling times for the accelerated leptons will be too high to maintain the power-law along 
the jets unless there is continuous distributed acceleration. Distributed, continuou s in sit u injection is a standard requirement for 
explaining the persistence of particle distributions along AGN jets (Je ster et al.l l2001) and seems to be a common feature in 
outflows. 

So to summarize the particle distribution discussion, we assume that the particles enter the jet with a relativistic quasi- 
Maxwellian distribution, the temperature of which, T e , is a fitted parameter. Some fraction of the particles are then assumed 
to be accelerated continuously beyond a location in the jets, z acc . Where exactly can be constrained by the frequency of the 
optically thick-to-thin break in the synchrotron spectrum, if detected, as well as from the synchrotron component contributing to 
the soft X-rays. 

Jet Continuum Components 

As the leptons travel outwards along the jets, they interact with the local magnetic field and photon fields and cool via syn- 
chrotron and inverse Compton radiation. Inverse Compton processes are strongest near the base of the jets, where the density 
is highest, and thus one X-ray component is due to the upscattering of both jet synchrotron photons as well as thermal disk 
photons by the quasi-thermal leptons in the base. For the assumptions in this paper (the thin disk is recessed with inner radius r m 
calculated from fitted parameters and 7j n ), the disk thermal photon field is only an important factor at the very base of the jets. 
This is because even the mild beaming in the jets quickly serves to weaken the disk photon field's energy density compared to 
the locally produced synchrotron photons. 

Once the acceleration zone is reac hed, it is difficult to suppress X-ray synchrotron radiation from the tail of accelerated leptons 
in the diffusive acceleration case (see lMarkoff et alJ2001al) . As described above, we are relaxing the constraint on the synchrotron 
cutoff by allowing for an alternate acceleration process. In this scenario, a power-law of synchrotron radiation will contribute 
significantly at least to the soft X-ray band, but not necessarily to the hard X-rays which are in this model dominated by inverse 
Comptonized disk and synchrotron photons. 

This distinction is important because, as shown by Markoff & Nowakl d2004h . a very typical hard state model where the X- 
rays are exclusively due to synchrotron radiation cannot easily reproduce disk reflection fractions greater than a few percent 
(assuming perpendicular geometry and no disk flaring). The larger distance of the first synchrotron emission region z acc from the 
disk, combined with moderate beaming, reduces its importance for disk feedback. On the other hand, inverse Compton processes 
near the base of the jets can produce reflection fractions of up to ~ 15-20%. This value is actually a lower limit because it 
does not in clude the effe cts of light bending, which will serve to increase the reflection fraction possibly quite significantly 
( iMiniutti & Fabianl2004l) 

Beyond relaxing the assumption about the acceler ation rate, the re a re a few other differences between this model and the 
synchrotron-dominated models considered in Markoff et al] ll2001al) and[Markoff et al. (2003). The comparatively weak inverse 
Compton component in prior models followed mainly from our choice of a small jet nozzle radius, ro ~ 3r g , which was based 
on the assumption that the jet radius was on the order of the event horizon. If, on the other hand, the jet base is contiguous with, 
or generated in, an extended corona, a larger scale seems more sensible. With this in mind, we here consider models with with 
typical values of ro > 10r g . 

The dependence of the calculated spectrum upon the model parameters, and their interdependence, are complex. Increasing the 
scale of the jet base decreases the electron density as well as the magnetic field, for a fixed equipartition relationship. This allows 
one to consider electron temperatures several times higher than those used in our previous models (to make up for lost synchrotron 
flux). The higher electron temperatures lead to greater inverse Compton emission relative to synchrotron processes in the X-ray 
band. Alternately, one can compensate for the weakened radiating power of the decreased particle density by increasing the total 
power input into the jet, Qj . Again, for a fixed equipartition relationship, the particle density is more sensitive to this change than 
the magnetic field (because the pressure of the field is oc B 2 ) and thus the inverse Compton component will experience a greater 
boost relatively. 
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Table A2. GX 339-4 and Cyg X-l Jet Model Fits (90% confidence level 

ERROR BARS) 



ObsID 


N s 

(10- 3 L Edd ) 


CIO" 3 L Edd ) 


ro 

(GM/c 1 ) 


(10'°K) 


p 


k 


p/f 






20181-01-02 
40108-02-01 
40108-02-03 


i 97+0.03 
L " L '-O.02 

64+ - 02 

u - DH -0.03 
34+0 01 
u - J +-0.03 


7.9 
4.7 
2 J 


20.2+0- 3 
9 6+ 05 
10.9+0;| 


5 01+ 07 
J -"'-0.08 
r 9O+0.13 
J - -0.12 

4.0 


9 Q4+0.01 
z -^-0.02 
9 9Q+0.01 

z, -''-o.oo 

2 65+0- 01 
-0.00 


1 7 ir+0.09 
1,,J -0.09 
1 19+0.01 

'■ lz, -o.oi 

2.2 


67+0- 03 

u - D '-0.02 
74+0-05 
u - /4 -0.01 

0.7 


25 
302 
87 


1 34+O.O2 
J -0.02 
1 41+O.OI 

1 99+0.03 


40099-01-19 
60090-01-26 
60090-01-41 


n oc+0.02 

74+0-0' 
v - '^-0.01 
n 70+0.02 


7.4 
7.1 
6.4 


o !+o.2 

4 4+0.2 

g 7+0.1 
o, '-0.3 


3 62+ 05 

T 9SJ+0.01 

3 - z!5 -o.oi 
3 90+° 03 

J - w -0.02 


2 50+ 003 
z,JU -0.02 

2 61+ 00 ' 
z - D1 -o.oi 

2 65+° 00 

ZOJ -0.01 


1 99+0.18 

1 77+0.04 

'■"-o.oi 
1 70+ 00 ' 

1 -' u -0.13 
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u, '°-0.O4 
73+O O8 
u ' /;> -0.02 

n c 1+0.03 
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16 
9 
/4 
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1- -"-0.02 
1 10+O.OO 
^^-O.OO 

1 17+0.05 

1-1 '-0.00 



NOTE. — We failed to resolve 90% confidence level error bars for parameters listed in italics. Jet model parameters are described in the paper and Appendix. Upper half of the table 
refers to GX 339—4, while the l ower half refers to Cyg X-l. We fixed t he mass, distance and inclination of GX 339—4 and Cyg X-l to 7 and 10 Mq, 6 and 2.5 kpc, and 30° and 47°, 
respectively i Hynes et al. 2003, 2004; Wu et al. 2001 ; Nowak et al. 1999). "Total output power from the jets in their rest frames L } is listed for completeness but is not itself a fitted 
parameter (see discussion in the Appendix). 



Table A3. GX 339-4 and Cyg X-l Jet Model Fits, Continued 



Obs ID 


"acc/f 




L d isk 
(10- 3 L Edd ) 


(keV) 


% 
(Tt) 


(io- 2 ) 


^line 

(keV) 


a 
(keV) 


Q/2tt 


X 2 /DoF 
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z -^-0.6 
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0.88+O;J0 


0.5 
0.1 
0.2 


4 5+0.4 
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z -'-0.5 
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02+ ' 
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NOTE. — We failed to resolve 90% confidence level error bars for parameters listed in italics. Jet model parameters are described in the paper and Appendices. Line parameters 
correspond to the usual XSPEC/ISIS normalization, centroid energy, and width of the gauss line model. Q/2ir is the reflection fraction from the XSPEC/ISIS reflect model. 
Upper half of the table refers to GX 339-4, while the lower half refers to Cyg X-l. ''Values listed for inner disk radius r m are derived from the spectrally difficult-to-constrain fitted 
parameters and Tdisk, and are thus also unconstrained: see discussion in Section|4| 



Table A4. GX 339-4 and Cyg X-l Coronal Model Fits (90% Confidence 
Level Error Bars) 



Obs ID 

(l(r 3 ) (10~ 5 ) (keV) (10~ 2 ) (keV) (keV) (4?rF/n) 

20181-01-02 0.21+O00 5.85+0;^ 2.86^; 2 5 4.0+^ 1.30+°;^ 0.14+0;^ 6.4+0;j 0.5+°:} 0.11+°^ 680+]«° 170/111 

40108-02-01 0.11+0;0 3 7.38+gjf 179+£f 1.35+0'f 0.06$$ 6.4+|];| 0.4+°; 3 0.10+°;^ 560+ 3 ™ 117/90 

40108-02-03 0.03+000 3.73^;^ 6.69^| 2 1.0+°;| 1.10+° 3 ° 0.02+°;<» 6.4+°; 2 0.7+°; 3 0.00+°; 2 ;> ... 17/37 

40099-01-19 9.23+^-JJJ 3.60+0;°' 1.35+°o> 422+f 0.85+°;°! 2~5+°] 6.2+0;' 08 +o.o o.l9+°-°| 1 140+}°° 226/158 

60090-01-26 6.78+0;|4 5.88+«;!]f 1.64+0; 02 47 + ' 2 0.87+°;' 3 1.9+°^ 6.0$J 0.9+°;° - 22 -o:o! 2+ -2 191/181 

60090-01-41 6.09+j;^ 5.67+0^ l-49+0;0 2 49+|J 0.85;0; 5 1.4+Oj 6.2+0;° .8+0;' 0.22+0;0j 2+' 5 186/183 



NOTE. — Model parameters are for the XSPEC/ISIS implementations of the eqpair model with ionized, smeared reflection, with additional contributions from a diskpn and 
gauss model. Upper half of the table refers to GX 339-4, while the lower half refers to Cyg X-l. 



